Ground state properties of a dilute homogeneous Bose gas of hard disks in two 

dimensions 
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The energy and structure of a dilute hard-disks Bose gas are studied in the framework of a 
variational many-body approach based on a Jastrow correlated ground state wave function. The 
asymptotic behaviors of the radial distribution function and the one-body density matrix are an- 
alyzed after solving the Euler equation obtained by a free minimization of the hypernetted chain 
energy functional. Our results show important deviations from those of the available low density 
expansions, already at gas parameter values x ~ 0.001. The condensate fraction in 2D is also 
computed and found generally lower than the 3D one at the same x. 



I. INTRODUCTION 

The study of two dimensional (2D) quantum many- 
body systems is a subject of considerable interest both 
from the theoretical and the experimental points of view. 
Atomic helium monolayers have attracted much atten- 
tion in these respects Q, and recently Bose-Einstein con- 
densation (BEC) in a strongly asymmetrical, quasi two- 
dimensional trap has been achieved^, 0- In a homo- 
geneous 2D gas the one-body density matrix acquires 
a power law decay below some low critical temperature, 
T < T c , 4] and, at T—0, its r — > oo limit differs from zero. 
This behavior defines a true condensate and BEC can 
occur in such a limit. The modified density of states in 
confined geometry, as it is the case for atoms in harmonic 
traps, makes BEC appear even at finite temperature in 
2D|1I1. 

In a recent paperQ, referred as I hereafter, we have 
examined the structure of the ground state of a 3D ho- 
mogeneous gas of bosons, interacting through both hard 
and soft core potentials. The study was carried out to 
enlighten the role of the interaction induced correlations 
along the density of the system and the possible occur- 
rence of universality in the dependence of different prop- 
erties on the gas parameter, x — na 3 (n being the particle 
density and a the s-wave scattering length). The calcu- 
lations were performed within the Correlated Basis Func- 
tions (CBF) theory, using optimized Jastrow correlated 
wave functions ||- The results were in excellent agree- 
ment with the existing Diffusion Monte Carlo[^ ones, ob- 
tained by the exact stochastic solution of the many-body 
Schrodinger equation. Moreover, the reliability (and lim- 
itations) of the analytical expansions in x of the energy 
per particle and of the condensate fraction were clearly 
assessed [To|. 



In this work we use the same variational approach to 
study the T = ground state of a homogeneous gas of 
N — ► oo hard-disks (HD) of radius a. We pay partic- 
ular attention to the long-range structure of the wave 
function, which is in turn intimately related to the long- 
range order of the density matrix and to the condensate 
fraction. 

At present, there are not published DMC results for 
these systems in the literature. However, based on the 
3D case experience we believe that the variational ap- 
proach may provide a very accurate description of the 
2D ground state properties and serve as a test of valid- 
ity for the available low density expansions (LDE). This 
fact may be of relevance also to the study of BEC in 
2D harmonic traps, since LDE are related to the Gross- 
Pitaevskii equation and its modification^! • LDE for 
the energy per particle in terms of the 2D gas parameter, 
x = na 2 , have been derived by several authors, starting 
from the leading orderfl2j: 
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(1) 



where the energy is in units of h 2 /2ma 2 . 

An interesting question, that we do not address in this 
paper, is the insurgence and the degree of universality 
in 2D. In the 3D Bose gas a dependence in the energy 
on the shape of the potential appears already at x ~ 
0.0010, El HH, and breaks down earlier for quantities 
other than the energy, as the short range structure of 
the distribution function and the condensate fraction. A 
similar study for the 2D gas would be of great interest by 
itself and also in view of the analysis of the quasi 2D BEC 
experiments in terms of the Gross-Pitaevskii equation. 

In the case of strong interactions, as a hard-core poten- 
tial, correlation effects are crucial. Within CBF theory 
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they are included by means of a many-body correlation 
operator, J-(l, 2, . . . , N), acting on the non-interacting 
ground state wave function, 4> (1, 2, . . . , V),||: 



*o(l,2, 



,iV)=^(l,2,...,iV)$ (l,2,...,iV), (2) 



with <&o = l for homogeneous bosons. T may be deter- 
mined according to the variational principle by minimiz- 
ing the ground state energy. 

In I a Jastrow correlation operator was considered, 



F=*j(l,2,...,N) = J} /(»"«): 



(3) 



where the two-body correlation function, /(ry), depends 
on the interparticle distance, r%j = \vi —Tj\, vanishes for 
m < a and goes to unity at large distances. Expecta- 
tion values of operators on the Jastrow correlated wave 
function correspond to multidimensional integrals (we 
consider only operators dependent on the spatial coor- 
dinates) which can be either directly evaluated by Monte 
Carlo type integrations (variational Monte Carlo, VMC) 
or expanded in cluster diagrams. The hypernetted chain 
(HNC) set of integral equations allows for summing an 
infinite number of diagrams, but it can be only approx- 
imately solved, since the class of the elementary dia- 
grams is not summable in a closed wav|l4j. The approxi- 
mation amounting to disregard the elementary diagrams 
(HNC/0) is expected to be reliable in the low density 
regime, as shown in I. 

In this paper we consider a system of N spinless bosons 
of mass 77i on a surface S, described by the hamiltonian 



h 

2m 



2 N 



N 



" = E v l + E ( 4 ) 

j=l l=i<j 



where V(r) is a two-body, symmetric hard-disk poten- 
tial, 



V(r) 



oo r < a 
r > a 



(5) 



and the radius a corresponds, for this potential, to the 
2D scattering length. We consider the system in the ther- 
modynamic limit (N and S — > oo, keeping the density, 
n = N/S, constant). 

Minimization of the energy with respect to f(r) pro- 
vides the optimal Jastrow factor, which can be ob- 
tained throu gh t he solution of the Euler-Lagrange (EL) 
equation[H ^3, 6E[f]/6f = 0. We adopt here the 
HNC/0 energy functional to solve the EL equation. The 
resulting correlation, fEL{r), shows a very long range 
structure (/el(t — > oo) — > 1 — a/r, with a constant) not 
easily accessible to VMC, due to the limited size of the 
simulation box. As a consequence we have also used a 
parametrized shorter range correlation factor, suitable to 
be used in VMC and whose parameters are variationally 
fixed. The comparison between the HNC/0 and VMC 
calculations with this short range correlation provides a 



check of the accuracy of the truncation in the cluster 
expansion. Diffusion Monte Carlo is also likely to suf- 
fer from the long range behavior of the wave function, 
as these fine structures, related to collective effects, are 
hardly distinguishable by numerical, finite size simula- 
tions. 

The plan of the paper is as follows: the correlated vari- 
ational theory and its implementation (HNC and Euler 
equations) are briefly outlined in Section II; an analysis 
of the EL asymptotic behaviors is presented in Section 
III; Section IV presents the results for the HD model; 
while summary and conclusions are given in Section V. 



II. VARIATIONAL THEORY 

Given the Jastrow correlated wave function ®, the 
energy per particle of the 2D homogeneous Bose gas is: 



E 1 

N = 2 n 



dr 12 g(n 2 ) 



V(r 12 ) 



2m 



V 2 ln/ 2 (n 2 ) 



(6) 

This expression simplifies for the HD potential J^J, 
since the radial distribution function (RDF), 



N(N-l) J dr 3 dr i ...dr N \^ \ 2 
n 2 J dridr 2 . ■ . drjv| v I'o| 2 ' 



(7) 



vanishes inside the core of the potential and only the 
kinetic part contributes to the energy. 

The HNC equations provide a procedure to evaluate 
9(r), 



N(r 12 ) 



fl(r 12 )e N ^ +E ^ 

77 / dr 3 [g(r 13 ) - 1] [5(7-32) 



1- N(r 32 )}(8) 



where N(r) and E(r) are functions representing the sum 
of the nodal and the elementary diagrams, respectively 
[l4|. The function E{r) is an input to the HNC equa- 
tions, which can be solved once a choice for it has been 
done. The simplest possible approximation corresponds 
to set E(r) = (HNC/0). This apparently drastic trun- 
cation is, however, justified at low densities since the ele- 
mentary diagrams, due to their high connectivity, do not 
appreciably contribute at the low densities relevant to 
BEC experiments. Otherwise, the energy and the RDF 
can be stochastically evaluated by Monte Carlo sampling 
of the corresponding many-body integrals JJJ. The ex- 
plicit calculation can be p erformed by using the standard 
Metropolis algorithm [T3 (see Ref. [Tg for a detailed de- 
scription of Monte Carlo methods). 

Connected to the RDF is the static structure function 
(SSF), S(k), 



S(k) = 1 + 77 / dre tk r (g(r) - 1) 



(9) 



often used in the analysis of (and extractable, within 
some approximations, from) scattering experiments in 
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condensed matter physics (e.g. neutron scattering off 
liquid Helium[l||). 

The Bose condensate is linked to the non-zero, long- 
range order of the one-body density matrix (OBDM), 



m(ri ,ri/) = ni(rn') 

J dr 2 ---dr J v*o(l,2 1 . 



N- 



,JV)*o(l',2,. 



(10) 

.,N) 



Jdr 1 ---dr N \^ Q \ 2 

In fact, Jii(rn/ — > oo)/n = hq, where no is the con- 
densate fraction. The depiction of no with respect to 
unity is an unmistakable indication of interparticle in- 
teractions (and, as a consequence, of correlations). The 
Fourier transform of the OBDM provides the momentum 
distribution (MD), 

n(k) = J dre ik,r m(r) (11) 
= (27r) 2 n n S(k) + / dr e ikr (m (r) - n no) 



As for the distribution function, the correlated OBDM 
can be computed by using HNC theory (2jj. In fact, n\(r) 
can be expressed in terms of new nodal and elementary 
functions as 



ni(r) 



n e 



N mw (r)+E ww (r) 



(12) 



where N ww (r) is the solution of a generalized HNC equa- 
tion. Again, setting E ww {r) = these equations can be 
solved in the HNC/0 approximation. 

An appropriate choice of the correlation factor is es- 
sential for the effectiveness of the variational approach. 
As stated in the Introduction, the best choice is the one 
satisfying the Euler-Lagrange equation, 



5E\g\ 
Sg{r) 



= 0. 



(13) 



which has been written in terms of the RDF rather than 
f(r), since in the HNC/0 scheme there is a one-to-one 
correspondence between these two quantities. The cor- 
relation function can then be obtained by inversion of 
the HNC equations. The EL equations are solved both 
in configuration and momentum space, as discussed at 
length in I for the 3D hard-sphere case. The theory can 
be straightfordwardly applied to the 2D hard-disks gas 
and it is briefly outlined in the Appendix. 

The correlation function produced by the solution of 
the EL equation shows a long-range structure that is 
discussed in the next Section. Finite size Monte Carlo 
techniques have difficulties to correctly deal with this 
long range behavior and correlations healing to unity in- 
side the simulation box are used. In this respect, we 
have also adopted a parametrized shorter-range corre- 
lation, /s_R(r), obtained by minimizing the energy com- 
puted at the two-body order of the cluster expansion, 



<?2_b(?~) = / 2 (0i constrained by a normalization condi- 
tion 



+/m / dr(l-/ 2 (r)) 



(14) 



which requires the use of a Lagrange multiplier /i. The 
minimization is performed under the healing condition 
fsR{r > d) — 1, while fi is fixed by imposing f'(d) = 
f"(d) = 0. The healing distance d is taken as a varia- 
tional parameter to minimize the many-body energy. 

For the particular case of the hard-core potential, the 
short range correlation, /s.fj(r), vanishes at r < a, while 

, ( < < n Y (Xa)Jo(Xr) - J (Xa)Y (Xr) 

fsR[a <r <d) = ^tttj; > U 5 ) 

Y (Xa)J {Xd) - J (Xa)Y {Xd) 

where A = —Am/i/h 2 and Jo and Y"o are Bessel functions 
of the first and second kind, respectively. In the d — > oo 
limit, fsR coincides with the exact solution of the 2D 
zero energy Schrodinger equation, f(r) tx ln(r) pi) . 



III. ASYMPTOTIC BEHAVIORS 

Applying to the 2D gas the sum-rule analysis of Rcf . 8] , 
it can be shown that the SSF at low momenta has a linear 
dependence, similar to the 3D gas: 



S(&->0) 



2mc ' 



(16) 



where c is the sound velocity in the medium. In terms of 
the RDF and of the correlation function, it corresponds 
to a long range behavior, different, however, from the 3D 
case where {g^D — 1) oc l/r~ 4 and (fan — 1) oc l/r~ 2 at 
large r-values. 

In fact, the 2D Fourier transform for a function, h(r), 
with circular symmetry can be written as: 



h(k) = lixn I drrf(r)Jo(kr), 



while h(r) is obtained from h(k) as 



1 



h(r) = — / dkkh(k)J (kr) 



(17) 



(18) 



Assuming that H(k) = kh(k) is a well behaved func- 
tion at the origin, one finds 



h(r — > oo) 



H(0) IH U (0) 3H IV {0) 



2 r 3 



+ ■ 



+ •••, (19) 



where only the k = values of the function and of its 
even derivatives enter in the expansion. 
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The RDF is obtained from S(k) through: 
1 



? (r) - 1 



2ttti 



dk k (S(k) - l)Jo(kr) , (20) 



so one readily obtains, from Eqs. (|16[) and (|19|) . the 
asymptotic limit 



g(r -> oo) = 1 



1 ft 1 

27m 2mc r 3 



(21) 



By inverting the HNC/0 equations (JSJ one finds for the 
nodal and correlation functions the following limits: 



N(r -> oo) 



1 2mc 1 
27m ft r ' 



and 



/(r ->• oo) -> 1 - 



1 2mc 1 

4/7m ft r 



(22) 



(23) 



showing that in 2D, correlations have longer range than 
in 3D. 

The long range structure of the OBDM is derived from 
the previous expressions and the HNC equations. In fact, 
given the structure we obtain for N^(r) 



N ww (r — > oo) 



1 mc 1 

Aim ft r 



and for the OBDM 

ni(r —>■ oo) 



n + n 



1 mc 1 
Aim ft r 



(24) 



(25) 



The momentum distribution has the same long- 
wavelength limit shown in 3D , namely 



lim kn(k) = — — 
fe^o w 2ft 



IV. RESULTS 



(26) 



In this section we present and analyze results for the 
energy, radial distribution function, static structure func- 
tion and one-body density matrix of the hard-disks gas. 
We have used the optimized Jastrow wave function ob- 
tained from the solution of the Euler-Lagrange equations 
and the short-range correlation of Eq. (|15fl , mainly to es- 
tablish a comparison between the HNC/0 and the VMC 
approaches. In the following, dimensionless quantities 
will be used: energies and distances will be given in units 
of ft 2 /2ma 2 and a respectively. 

Several corrections to Elo have been proposed in the 
literature. Kolomeisky and Straley [2j| used renormaliza- 
tion group techniques to study the ground state of dilute 
Bose systems as a function of the space dimensionality. 
They found a general expression, valid for strong inter- 
actions when x — > 0, that simplifies for hard disks to 



E 



KS 



Elo 



I Ms) I 

ln(47rx) 



1 - 



In (- ln(47rx)) 



ln(47rx) 



(27) 



Cherny and Shanenko \2A 
pansion, 



derived an alternative ex- 



E, 



(u) 

cs 



E 



LO 



u 



in the parameter u satisfying the equations: 



u = (5(1 + it In u) 



5 = 



1 



ln(7rx) + 27 ' 



(28) 



(29) 



where 7 = 0.577 ... is the Euler constant. These authors 
also gave an expansion of u in terms of 5, allowing to 
write the series (|28f) in the form: 



E. 



cs 



E 



LO 



I ln(x) 



S + 5 2 In 6 + 



S l+S 3 \n 2 6 



25 A In 5 



Table [I] reports the energy per particle as a function 
of x in the EL approximation (EL/HNC), the variational 
Monte Carlo calculation starting from fsn(r) (SR/VMC) 
and the HNC approach with the same correlation func- 
tion (SR/HNC). We remind that the HNC/0 approxima- 
tion is used everywhere but in the VMC. The results of 
the x-expansions previously discussed are also reported. 

The comparison between SR/VMC and SR/HNC 
shows that the influence of the missing elementary di- 
agrams on the energy is less than 1%, except at the 
highest value x = 0.1 (~ 2.3%). This gives us confi- 
dence that the variational principle is mostly satisfied 
within our HNC/0 calculations, providing the hierarchy: 
Eexact < Eel < Esr. Only at x = 0.1 these inequal- 
ities do not numerically hold. The cases x = 10~ 5 and 
x = 5 • 10~ 2 can be considered to fulfill the inequali- 
ties if we take into account the numerical accuracy as- 
sociated to the calculation at these quantities. If we es- 
timate the contribution of the elementary diagrams on 
Eel by scaling it by the ratio E^Jf c /Eg£ c , we obtain 
Eel/N(x — 0.1) = 0.9075, restoring all the inequalities. 

Lieb [25j pointed out that a lower bound to the exact 
energy is given by 



Eiow — E 



LO 



l-O 



(l ln(x) 



-1/5 



(31) 



and that E exact / Elo ~ * 1 when x — > 0. Both the varia- 
tional energies (EL and SR) comply with condition l|31|) 
and seem to tend to E exact when x goes to zero. 

Table JIJ gives also the healing distance, d, of the SR 
correlation in units of a. d increases when x — + 0, the La- 
grange multiplier decreases and the energy goes to zero. 
Therefore, fsR(r) x ~>o can be approximated by its A = 
limit, which coincides with the zero energy limit of the 
two-body Scrodinger equation, 



Mr) 
ln(d) 



(32) 



These results are also shown in Figure QJ, where the 
variational scaled energies per particle (EL, SR/VMC 
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and SR/HNC) are compared with the Eks and Ecs es- 
timates. All energies have been divided by E^q in order 
to stress the deviations from the low-density limit. The 
limit is approached by Eel and E$r from above when 
x decreases, although it has not been yet fully reached 
at x = 1CP 5 . The differences between the variational 
and the low density energies are still visible^even in the 
density range relevant to BEC experiments |3j. 

Eks does not satisfy the x = lower bound and, 
starting from x ~ 0.005, becomes higher than the vari- 
ational upper bound provided by Eel and Esr- Ecs 
satisfies the low density limit, but lies above the varia- 
tional upper bounds at any value of x. ^ s always 

larger than Eq S and the difference increases drastically 
along x. Notice that Eq. H29[) does not have solution at 
x > 0.0369. 

The EL optimization procedure does not significantly 
affect the energies obtained with fsR^), since the en- 
ergy is dominated by the short range structure of the 
HD potential, requiring g(r) to vanish inside the core. 
The effects of the long range structure of the EL corre- 
lation are, instead, clearly evident in the behavior of the 
radial distribution function. 

The EL RDF is shown in Figure J5| for different values 
of x. At low x, <7Bi(r) is a monotonically increasing func- 
tion of the distance. However, it develops a local maxi- 
mum close to the core radius at densities x > 0.01. This 
is a genuine many-body effect induced by the strong cor- 
relations at high density. The same behavior was found 
in I for the 3D Bose gas. As expected, the correlation 
hole is more pronounced at larger densities. 

The long range limit of §el{t) is shown in Figure at 
x = 0.01 and x = 0.1. The quantity shown is ir 3 (j(r) — 
1) whose dimensionless asymptotic limit is: 



^ 2ttc ' 



(33) 



This ratio is smaller at x = 0.1, implying that the sound 
velocity increases with expected. Consistent with 

the previous figure, the asymptotic limit is reached faster 
at larger densities. Figure gives the EL correlation 
and the nodal function, N(r), at x = 0.001. We show 
[xr(l — /Ei(r))] and [xrN(r)] to enlighten the asymp- 
totic limits <|23|) and (|22|l whose dimensionless values are 
c/47r and c/2-7r, respectively. The fact that one limit is 
twice the other is clearly appreciated in the figure. Notice 
also that due to the chain process implied by the HNC 
scheme, N(r) merges to the 1/r law at larger distances 
then f(r). To illustrate the different asymptotic behav- 
iors, we also show [x r(l — fsR(f))] which goes quickly to 
zero. 

The EL static structure function, S'(fc), is shown in 
Figure At low densities the SSF reaches the asymp- 
totic value, S(k) — > 1, already at k ~ 1. As in the RDF 
case, the overshooting of the SSF at the highest density, 
x = 0.1, is a consequence of the correlations. The linear 
regime of S{k) around the origin is appreciable, although 



the calculation of the ratio S{k)/k shows deviations from 
a constant value already at low k. 

In figure © we plot the one-body density matrix, 
rti(r), in the EL approach for x = 0.01, 0.005 and 
0.001. The asymptotic limit of ni(r) defines the value 
of the condensate fraction, which decreases when the gas 
parameter increases. The asymptotic value is reached 
faster when x increases. The detailed long-range behav- 
ior (|2*5|l is presented in Figure J7J by showing the quantity 
r [ni (r) / {n^n) — 1] , whose dimensionless asymptotic value 
is c/(8irx). Even if the speed of sound increases with 
x, the value of this limit is dominated by the presence 
of the gas parameter in the denominator and the over- 
all quantity increases when x decreases. Finally, the EL 
and SR/VMC condensate fractions, uq(x), are reported 
in Figure ©■ The Figure also contains the low-density 
prediction [12J, 



LD 



ln(a;) 



(34) 



,LD 



appears to sensibly overestimate the condensate 
fraction. The EL and SR/VMC condensate fractions are 
very similar except for the largest value of x reported in 
the figure, where the contribution of the elementary di- 
agrams could be important. This fact indicates that the 
value of no is not very much affected by the inclusion of a 
long-range structure into the correlation function. How- 
ever, the use of /bl(^) is crucial to approach this value in 
a proper way, that is, to satisfy Eq. lj2*5l) . Also reported in 
the figure is the condensate fraction of the 3D system of 
hard spheres, taken from I. At fixed x, the 2D condensate 
fraction is smaller than the 3D one, indicating that cor- 
relations in the 2D system are stronger, thus promoting 
more particles outside the zero-momentum state. 



V. SUMMARY AND CONCLUSIONS 

In this work we have analyzed the energy and structure 
of a homogeneous gas of bosons in two dimensions inter- 
acting via a hard disk potential whose core radius equals 
its corresponding 2D scattering length. We have adopted 
a variational many-body approach, based on a Jastrow 
correlated ground state wave function. The expectation 
values have been computed both in the framework of the 
hypernetted chain theory (within the HNC/0 approxi- 
mation) and with the variational Monte Carlo method. 
Two types of correlation functions have been used: (i) a 
long range one, obtained by the free minimization of the 
HNC/0 ground state energy and preserving the correct 
asymptotic behaviors of the wave function; (ii) a short 
range one, to be used in the Monte Carlo sampling and 
providing a check of the accuracy of the cluster expan- 
sion. 

By comparing with the VMC results, the accuracy of 
the HNC/0 energies are better than 1%, except at the 
highest density, x — 0.1, where the error is still less than 
3%. The EL minimization lowers the energy with respect 
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to the SR correlation by ~ 1.5%. We do not expect fur- 
ther large reductions from a complete DMC calculation. 
The low density expansions start to severely deviate from 
the variational results already at x ~ 0.001, and the most 
accurate of them appears to be the Cherny and Shancnko 
expansion in terms of the parameter 5. However, their 
use for estimating corrections to the 2D Gross-Pitaevskii 
equation, especially in the large gas parameter regime, 
seems questionable. 

Finally, the condensate fractions lies well below the val- 
ues predicted by the low density theories and also below 
the results for the three dimensional gas of hard spheres 
at the same gas parameter. 

We conclude that the variational theory is a power- 
ful and reliable tool to study dilute systems, also in 2D. 
Moreover, the homogeneous gas results may be used in a 
local density type approximation |Tlj to analyze bosons 
in two dimensional harmonic traps. 



VI. APPENDIX 

In this appendix we discuss the EL equations for an 
interacting many-body system in 2D. As in the 3D case, 
the solution to the optimization equation, 



5E\g\ 
Sg{r) 



= 0. 



can be obtained in momentum space, yielding 



S(k) 



t(k) 



(35) 



(36) 



where t(k) = h 2 k 2 /2m is the free-particle energy spec- 
trum and V p h{k) is the particle-hole interaction. The 
latter can be expressed in configuration space and reads, 
disregarding the contribution of elementary diagramas 



V ph (r) = g(r)V(r) + — V^g(r) + [g(r) - 1] wj(r) , 

(37) 

in terms of the induced interaction cj/(r). In momentum 
space, we have 



1 „J2S(k) 



l)(S(fc)-l) 



S(k) 



(38) 



Eqs. H36fl to l|38|) are to be solved simultaneously. This 
can be done starting from a suitable choice for g(r), per- 
forming its FT to get S(k), evaluating u>i(k) and V p h(r), 
and then deriving a new S(k) with the help of Eq. 1)36(1. 
This procedure is iterated until the difference between 
two consecutive iterations is as small as desired. 

Up to this point there are no formal differences between 
the 2D and 3D cases. The main deviation is the way in 



which the Fourier transforms are carried out. In 2D and 
for a general function, /(r), the FT and its inverse read 



f(k) = 27m / drrf(r)J (kr) 



(39) 



f(r) = ^ J dkkf(k)J Q (kr) (40) 

where n is the (constant) density and Jo{x) the zero or- 
der Bessel function of the first kind. One way to im- 
plement these transformations is to use a finite box in 
configuration and momentum spaces of length L and K, 
respectively. The additional conditions f(r = L) = and 
f(k = K) = lead to a discretized set of coordinates and 
momenta, related to the zeros Xj of Jq(x) through the re- 
lations 



hj — —f~ i T a — L — 



JV 



(41) 



with j,a = 1, 2, . . . , N, N being the total number of 
points in the grids. A Gauss integration rule based on se- 
ries expansion in Bessel functions and the orthogonality 
relation, 



drrJ (kir)J (kjr) 



can then be built, leading to 



h > (42) 



/ drrf(r)J (kjr) = ^ u a Ja( k j r a)f(r a ) , 

/ dkkf(k)J (kr a ) = y2"jMkjr a )f(kj) , (43) 
,=i 



whith the integration weights 
2L 2 



uJ t , 



(44) 



Eqs. (|43|l turn integrals into algebraic products that 
can be carried out numerically in a neat and fast way 
using available standard libraries. 
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TABLE I: Energy per particle for the hard disk model, as a function of x. Eel is obtained by solving the HNC/0 EL equation, 
Evmc and Esr are the VMC and HNC/0 energies with the SR correlation having a healing distance d, in units of the scattering 
length, evmc is the statistical error in the variational Monte Carlo calculation. The last four columns give the energies in 
different low density expansions (see text). 
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FIG. 1: Scaled energy per particle of the HD gas as a function 
of x. Solid circles: EL/HNC; open squares: SR/HNC; open 
diamonds: SR/VMC; stars, pluses and crosses: low density 
expansions, KS, CS (u) and CS (S) , respectively. Notice that 
diamonds are hardly distinguished from the solid circles. 
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FIG. 2: EL radial distribution functions g(r) for 2D hard disks 
at several values of x. Solid, dot-dashed, long-dashed and 
short-dashed lines stand for x = 0.1,0.01,0.001 and 0.0001, 
respectively. 
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FIG. 5: EL static structure factor, S(k), at several values of 
x. Solid line: x = 0.1, dashed line: x = 0.01, dot-dashed line: 
x = 0.001, and dotted line: x = 0.0001. 




FIG. 6: EL one-body density matrices at several values of x. 
Solid line: x = 0.01, dashed line: x = 0.005, and dotted line: 
x = 0.001. Open circles, SR/VMC results at x = 0.005. 
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FIG. 7: Long range behavior of the EL one-body density 
matrix at several values of x. Solid line: x = 0.01, dashed 
line: x = 0.005, and dotted line: x = 0.001. 
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FIG. 8: Condensate fraction as a function of x. Black circles, 
triangles and solid line correspond to the EL/HNC, SR/VMC 
and low-density expansion results, respectively. Open squares 
stand for the corresponding 3D values taken from I. 



